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Abstract 



Since the thermal history of the Universe is unknown before the epoch of primordial nucleosyn- 
thesis, the largest temperature of the radiation dominated phase (the reheating temperature) might 
have been as low as 1 MeV. We perform a quantitative study of supersymmetric dark matter relic 
abundance in cosmological scenarios with low reheating temperature. We show that, for values of 
the reheating temperature smaller than about 30 GeV, the domains of the supergravity parameter 
space which are compatible with the hypothesis that dark matter is composed by neutralinos are 
largely enhanced. We also find a lower bound on the reheating temperature: if the latter is smaller 
than about 1 GeV neutralinos cannot be efficiently produced in the early Universe and then they 
are not able to explain the present amount of dark matter. 
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I. INTRODUCTION 



Supersymmetric dark matter provides one of the hottest topics at the border between 
Cosmology and Particle Physics. This is due to the fact that in i?-parity conserving Super- 
symmetric theories the Lightest Supersymmetric Particle (LSP) is stable and may provide 
the cold dark matter, whose existence is inferred by a large number of independent observa- 
tions [[1], 0]. Among the different supersymmetric candidates, the neutralino turns out to be 
a perfect dark matter particle, since it has neither charge nor colour, its only interactions 
being of the weak type. 

The present abundance of neutralinos depends on the thermal history of the Universe. 
In the early Universe interactions may keep neutralinos in thermal equilibrium with the 
radiation bath until their abundance freezes out at a temperature Tp. The neutralino mass 
is constrained by accelerator data to be heavier then a few tens of GeV. This implies that 
it decouples in the early Universe when it is non relativistic, at Tp in the GeV range. This 
picture is correct if the maximum temperature in the radiation-dominated era, which from 
now on we will refer to as the reheating temperature Trh, is much larger than the freeze-out 
temperature Tp. If this is the case, the neutralino relic abundance turns out to be: 

s P ^ k 2.6.10-°GeV-3 

Pc (tfann^) 

where p c = 1.8791/i 2 x 10~ 29 g cm -3 is the critical density, h is the Hubble constant in 
units of 100 km sec -1 pc -1 , while cr ann is the WIMP pair-annihilation cross section, v is 
the relative velocity between the two annihilating particles, and brackets denote thermal 
average. Indeed, what specifically makes the neutralino an ideal dark matter candidate is 
that in Eq. (j^) the value of the annihilation cross section, calculated in a wide variety of 
susy models, may yield a result for Q x h 2 which falls in the correct interval suggested by 
present day observations for the amount of non-baryonic dark matter in the Universe [|], [|: 

0.05 < n M h 2 < 0.3. (2) 

When exploring the neutralino parameter space and its chances of discovery both in ac- 
celerator and non-accelerator searches, this argument has been usually turned the other 
way around, and Eq. @ used as a constraint on the neutralino parameter space. Depend- 
ing on the particular supersymmetric scenario, this may have important consequences on 
the allowed supersymmetric configurations. In particular, in Super gravity-inspired models 
(SUGRA), the allowed neutralino parameter space turns out to be severely reduced by re- 
quiring that tt x h 2 falls inside the interval defined by Eq. (@) |, 1 1, |, ^ g |, 0, 0,|IJ,|lJ. 

The robustness of these constraints relies on the cosmological assumptions that lead to 
Eq. ([I]). Indeed, the thermal history of the Universe before the epoch of nucleosynthesis is 
unknown. The maximum temperature in the radiation-dominated era Trh may have been 
as low as 1 MeV (but not smaller in order not to spoil the nucleosynthesis predictions). The 
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possibility of a low reheating temperature of the Universe has been recently discussed in Ref. 
|Tj| . There it was shown that a low reheating temperature has important implications for 
many topics in cosmology such as axion physics, leptogenesis and nucleosynthesis constraints 
on decaying particles. In particular, it was shown that stable weakly interacting massive 
particles may be produced even if the reheating temperature is much smaller than the freeze- 
out temperature of the dark matter particles, T RH < T F , and that the dependence of the 
present abundance on the mass and the annihilation cross section of the dark matter particle 
may differ drastically from standard results 1 . 

The goal of this paper is twofold: first, we wish to extend the analysis of Ref. [1 I 



and perform a quantitative study of the case of neutralinos in SUGRA scenarios, analyzing 
in detail the impact that a low Trh may have for the present neutralino relic abundance; 
secondly, we aim at providing a lower bound on the reheating temperature. The logic is the 
following. All matter is produced at the end of inflation [[17]] when all the vacuum energy 
stored into the inflaton field is released and the Universe becomes radiation-dominated with 
the initial temperature T RH . During the reheating process, particles are generated through 
thermal scatterings and quickly thermalize. Among them, dark matter particles may be also 
produced but their final number depends strongly on the reheating temperature. If the latter 
is too small, the thermal bath does not give rise to a number of neutralinos large enough to 
make them good candidates for dark matter. This leads to a lower limit on Trh- We will 
find that the reheating temperature needs to be larger than about 1 GeV for neutralinos to 
be good dark matter candidates 2 . 

The plan of the paper is as follows. In Section II we briefly recall the calculation of a 
WIMP relic abundance both in the low-reheating scenario and in the standard case. In 
Section III the specific supersymmetric models which will be considered in our analysis are 
introduced. In Section IV we will discuss our results. Section V is devoted to our conclusions. 



II. NEUTRALINO RELIC ABUNDANCE: THE STANDARD SCENARIO AND 
THE LOW REHEATING TEMPERATURE SCENARIO 

In this section we outline the ingredients which are relevant to the calculation of the 
neutralino relic abundance both in the standard radiation-dominated scenario and in the 



low reheating early Universe scenario introduced in Ref. JLJ]. We address the reader to 
Refs. [|14|, |19|] for further details. 



1 Low reheating scenarios lead as well to a new perspective on baryogenesis 15 and to the possibility that 
massive neutrinos may play the role of warm dark matter pl| . 

2 In this paper we suppose that neutralinos are produced during the reheating process only through thermal 
scatterings. Another source might be the direct decay of the inflaton field into neutralinos Oj. This 
introduces though another unknown parameter, the decay rate of the inflaton field into neutralinos, and 
we do not consider this possibility any further. 
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A. The Standard Scenario 



The number density n x of neutralinos in the early Universe is governed by the Boltzmann 
equation which takes into account both the expansion of the Universe and the neutralino 
interactions in the primordial plasma: 

^ = -3Hn x - (a ann v) [n\ - (r^) 2 ] , (3) 
where H is the Hubble parameter, t denotes time and n eq is the neutralino equilibrium num- 

A 

ber density. In the r.h.s. of Eq. (|3]), the first term describes the Universe expansion, while 
the second term takes into account the change in the x number density due to annihilation 
and inverse-annihilation processes. 

At high temperatures, the evolution of n x closely tracks its equilibrium value n eq . In 
this regime the interaction rate of the \ particles is strong enough to keep them in thermal 
equilibrium with the plasma. As temperature decreases, for heavy particles like neutrali- 
nos the quantity n x q becomes exponentially suppressed and therefore the interaction rate 
T = n x (<T ann f) turns out to be rapidly ineffective in maintaining neutralinos in thermal 
equilibrium: when the x' s mean free path becomes of the order of the Hubble scale, x' s 
interactions are freezed-out and the x' s number density per comoving volume is freezed-in. 
This situation occurs at a temperature Tp (freeze-out temperature) and clearly depends on 
the strength of interactions. 

By integrating the Boltzmann equation of Eq. (D) up to the present time, one finds the 
neutralino relic density Q x h 2 : 

n x h 2 = m ^ {T = 0) h\ (4) 

Pc 

A simple analytic approximation of the solution of the Boltzmann equation ([J) allows us 
to write down explicitly both the neutralino relic abundance: 

Q x h 2 ~ 8.77 x KT 11 — ^ GeV " 2 . (5) 

gj (Tp }St d) (°ann^)int 

and the value of the freeze-out temperature for the standard scenario Tp = Tp tS td in the 
implicit form: 



xp.td ^ In 



0.038 * * F F f ((W); 

g* (Tp tSt d) 



(6) 



where xp^ st d = fn x /Tp^ std and (<J ann v)p denotes the value of (cr ann v) at the freeze-out tem- 
perature. In Eqs. (|5|||) Mp denotes the Planck mass, g is the number of internal degrees 
of freedom of x, 9*(Tp) is the effective number of degrees of freedom of the plasma at the 
freeze-out and (cr ann t> )i nt denotes the integrated value of ((T^v) from Tp up to the present 
temperature. Making use of the non-relativistic first-order expansion of (a ann f ) in terms of 
the variable x = m x /T: 

(cr^nnv) = a + bx' 1 , (7) 
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we have: (cr ann f)i nt = axp 1 + bx F 2 /2. Eq. (||) shows the well know result that the present 
abundance of a cold relic particle is inversely proportional to its annihilation cross section: 

A crucial point in this discussion is that, in the standard cosmological scenario, freeze-out 
occurs in a phase of the evolution of the Universe when the expansion is adiabatic and the 
energy density is dominated by radiation: T ~ a -1 and H ~ T 2 ~ a~ 2 . These relations 
between the temperature, the scale factor a and the Hubble parameter are modified in the 
reheating phase of the low reheating-temperature scenario discussed in the next Section: if 
freeze-out occurs during the reheating phase, a lower neutralino relic abundance at present 
time is obtained. 



B. A low reheating temperature scenario 



It is by now accepted that during the early epochs of the Universe there was a primordial 
stage of inflation |17| responsible for the observed homogeneity and isotropy of the present 
Universe as well as for the generation of the cosmological perturbations. 

The radiation-dominated era of the Universe is usually assumed to be originated by 
the decay of the coherent oscillations of a scalar field, the inflaton field, whose vacuum 
energy has driven inflation 3 . The decay of the scalar field into light degrees of freedom and 
their subsequent thermalization, called reheating, leaves the Universe at a temperature Trh-, 
which represents the largest temperature of the plasma during the subsequent radiation- 
dominated epoch, when temperature is a decreasing function of time. The onset of the 
radiation dominated era is in fact placed at the temperature Trh, at the end of the 
reheating phase. 

Usually Trh is assumed to be very large and - in any case - larger than the neutralino 
freeze-out temperature Tp. This fact implies that the present-day relic abundance of any 
particle which freezes-out at a temperature Tp < Trh is not affected by the history of the 
Universe during the reheating phase. However the only information we have on the smallest 
value of Trh is from requiring a successful period of primordial nucleosynthesis, T RH > 1 
MeV. Therefore, from a phenomenological point of view, Trh is actually a free parameter. 
This implies that the situation in which a relic particle decoupled from the plasma before 
reheating was completed (i.e.: T F > T RH ) could be a viable possibility, with important 
implications in the calculation of the cosmological abundance of relic particles. 



3 Identifying this scalar field with the inflaton field is not strictly necessary. It might well be identified 
with some massive nearly stable particle, such as some light modulus field present in supersymmetric 
and (super)string models. In such a case, after inflation the Universe might have been matter dominated 
by the energy density of this modulus and then become radiation dominated after its decay. In other 
words, there might have been more than one reheating process during the thermal history of the Universe. 
Needless to say, the one relevant for us is the latest. 
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Let us consider the scenario of early Universe discussed in Ref. ||14j| . During the reheating 
epoch, the energy density of the Universe is dominated by the coherent oscillations of a scalar 
field (p. This period begins at a time Hj 1 and lasts until a time T^ 1 set by the scalar field 
decay rate T^. The dynamics of the system for Hj 1 < t < T^ 1 is described by the Boltzmann 
equations for the energy densities p<p,R, x of the three coupled components: the (unstable) 
massive field 0, radiation R, and the (stable) massive WIMP x's: 



dt 
dp R 

dt 

dn x 

~dT 



SHptp - , (8) 
-4Hp R + I> + (2(E x ))(a ann v) [n\ - [n^f] , (9) 
-3Hn x - (<w) \nl - (n^) 2 ] , (10) 



where the quantity 2(E X ) represents the average energy released in each XX annihilation. 
Notice that here we assume that decays into radiation, but not into x's. in the following 
we will recall the main properties of the system described by Eqs. (§,[|,|10|), addressing the 
reader to Ref.jnj for a complete discussion. 

At very early times the energy density of the Universe is dominated by the scalar filed 
0, and both the radiation and WIMP energy densities are negligible. As the scalar field 
decays the temperature T grows until it reaches a maximum value Tmax and then decreases 
as T oc a -3 / 8 up to the temperature T RH at the time t ~ T^ 1 which determines the end 
of reheating (in this second stage the temperature of the radiation produced in the first 
stage is cooled down by expansion, and the entropy release due to the decay of induces 
a gentler cooling compared to the radiation-dominated case). This non-standard relation 
between the temperature and the scale factor may significantly affect the calculation of the 
X particles relic abundance, depending on the duration and on the details of the reheating 
phase. 

The relevant mass scales which can help in understanding the different regimes which 
may occur are: the mass of the WIMP m x , the temperature Tp of the WIMP freeze-out 
(which can happen before or after the reheating has been completed), the temperature at the 
end of the reheating phase T R h and the maximal temperature reached during the reheating 
phase Tmax- Two hierarchies are present for these mass scales: T R h < Tmax and, since we 
are dealing with cold relics which decouple when non-relativistic, Tp < m x . 

The relation between m x and Tmax determines whether the WIMPs, which are produced 
during the reheating phase, become relativistic {Tmax > rn x ) or not (Tmax < m x )- More 
important is to determine whether the x particles reach thermal equilibrium during the 
reheating phase. This is determined by the strength of their interactions, and in particular 
it depends on (cr ann t>). According to the values of the WIMP mass and annihilation cross 
section, two possible non-standard regimes may be schematically singled out for the WIMP 
relic abundance: 

(i) The x particles generated during reheating are always non-relativistic and never reach 
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thermal equilibrium. Integration of Eqs. (f^|j^|10"l) shows that the process of x production takes 
its main contribution around the temperature T* ~ m x/^ an d is exponentially suppressed 
outside a narrow interval centered on T* |I3 |. So most of the x particles are produced at 
T*. For T < T* the total number of x's is frozen and their density is diluted by expansion. 



The condition n x (T 
the order of 



0) < nf^X 1 *) implies an upper limit on (a ann t>)* = (c7 aim 'u)(T=T») of 



(aannV)* < 7 X 10~ i4 - 



10 



10 



1/2 



77t-i 



100 GeVy V T f 



/100 MeVV 



V 



GeV 



-2 



(11) 



In this regime the WIMP relic density turns out to be proportional to the WIMP self- 
annihilation cross section: Q x h 2 ~ (cr ann t> )*, instead of being inversely proportional, as in 
the standard case. We anticipate here that for the case of neutralinos in supersymmetric 
models, this situation only occurs when the reheating temperature is smaller than about 
300 MeV. In this situation the neutralino relic density is significantly suppressed. 

(ii) The x particles reach thermal equilibrium, and then freeze-out when non-relativistic 
before the reheating phase is concluded, i.e. at a temperature Tp = T F ^h > T RH . During 
the phase when they reach thermal equilibrium, the x particles may or may not become 
relativistic, depending on the value of Tmax- In any case, they decouple as non-relativistic, 
leading to a cold relic. The usual freeze-out condition is modified compared to the standard 
case because the energy density is dominated by the scalar field and the relation between the 
Hubble constant and the temperature is given by H oc T 4 , as compared to the radiation- 
dominated case where H oc T 2 . The first consequence of this fact is that freeze-out occurs 
earlier and the WIMP density at Tp^h turns out to be higher. However, as the Universe 
cools down from T Frh to T RH , due to the entropy produced by the decays, the expansion 
dilutes n x by a factor {T R H/T F}r h) s , which is much smaller than the dilution factor in the 
radiation-dominate case (Trh /Tp^td) 3 ■ When the two effects are combined, the final result 
is a suppression of Q x by roughly a factor T^ H T Fstd / (Tp )rh ) 4 as compared to the standard 
case. An analytic approximation of the ensuing result for the x re li c abundance is given by 



14 



Q^h 2 ~ 2.3 x 10 



-n9* 



1/2 (Trh) 



T 3 RH GeV- 2 



g*(TF t rh) m x (ax F 4 rh + 4bx F 5 rh /5) 



(12) 



where the freeze-out temperature during the reheat ing-phase T F ^ rh (x F>r h 

M 

/2 {T RH ) M P T RH : , 2 - A 2 , 

[axj h + 5bxj h /4) 



m x /Tj 



XF,rh 



111 



0.015 



1/2 

9 9* 



g*(Tj 



F,rh 



77K 



F,rh) IS 



(13) 



12Q and ([T3|) we have used the annihilation cross-section expansion of Eq. 



In both Eqs. 
(0)- 

Obviously, if the decoupling of the x particles from the plasma occurs after the reheating 
phase is concluded, the standard scenario is recovered and the relic abundance has the 
ordinary expression of Eq. (H) . 
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In the following Sections we will perform a detailed calculation of the neutralino relic 
abundance in a Supergravity framework in order to study the consequences of the low 
reheating temperature scenario outlined above. In our analysis we will make use of the 
analytical solutions given in Eqs.(|,||) and Eqs. (|T2l , |T3D , appropriately interpolated in the 
intermediate regime where Tp is close to Trh (the interpolation has been determined on 
the basis of numerical solutions of the relevant differential equations). We have numerically 
verified that the analytical solutions are accurate enough for our purposes when compared 
with full numerical solution of the relevant Boltzmann equations of Eq. ([|) and Eqs. (3,P,F0|) ■ 



III. THE NEUTRALINO IN MINIMAL SUPERGRAVITY 



Supersymmetric theories naturally predict the existence of viable dark matter candidates 
if .R-parity is conserved, since this symmetry prevents the lightest supersymmetric parti- 
cle (LSP) from decaying. The nature and the properties of the LSP depend on the way 
supersymmetry is broken. In models where supersymmetry breaking is realized through 
gravity- (or also anomaly-) mediated mechanisms, the LSP turns out to be quite naturally 
the neutralino, defined as the lowest-mass linear superposition of photino (7), zino (Z) and 
the two higgsino states (H±, H%): x = «i7 + (I2Z + a^Hl + a^H^. 

Even assuming a minimal supersymmetric extension of the Standard Model, the super- 
symmetric theories may be explored in a variety of different schemes, ranging from those 
based on universal §, § §, §, 0, g [[§ 0, 0, |TB[ or non-universal §, § §, §, |TT|, [T2| Super- 



gravity, where the relevant independent supersymmetric parameters are defined at a grand 
unification scale, to effective supersymmetric theories which are defined at the electroweak 
scale [§, |20], plfl . In the present paper we will mostly concentrate on the discussion of uni- 
versal Supergravity and we will comment on the results which can be obtained in different 
supersymmetric schemes. 

The essential elements of a generic minimal supersymmetric model are described by a 
Yang-Mills Lagrangian, by the superpotential, which contains all the Yukawa interactions 
between the standard and supersymmetric fields, and by the soft-breaking Lagrangian, 
which models the breaking of supersymmetry. Implementation of this model within a Super- 
gravity scheme leads quite naturally to a set of unification conditions at a grand unification 
scale (Mqut) for the parameters of the theory: 

• Unification of the gaugino masses: 

Mi(M GUT ) = TTli/2 , (14) 



• Universality of the scalar masses with a common mass denoted by mo: 

mi(M GUT ) = m , (15) 
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• Universality of the trilinear scalar couplings: 

A\M GUT ) = A d (M GUT ) = A U (M GUT ) = A m . (16) 

We denote this scheme as universal SUGRA (or minimal SUGRA). The relevant param- 
eters of the model at the electroweak scale are obtained from their corresponding values 
at the Mqut scale by running these down according to renormalization group equations. 
By requiring that the electroweak symmetry breaking is induced radiatively by the soft su- 
persymmetry breaking, one finally reduces the parameters of the model to five: mi/2, mo, 
Aq, tan/? ee 1)2/ V\ and sign(/i), where V\ and t>2 denote the vev's of the two Higgs field of 
the model and /i is a mixing parameter between the two Higgs fields which enters in the 
superpotential. These parameters are a priori undetermined. However, bounds coming from 
supersymmetry and Higgs searches at accelerators and on supersymmetric contributions to 
rare process, like the b — > s + 7 radiative decay, introduce limits on the model parameter 
space. Also theoretical arguments concerning the naturalness of the theory may be used in 
order to identify typical scales beyond which the main attractive features of supersymmetry 
fade away. For instance, fine-tuning arguments may be invoked to set bounds on m and 
7711/2 [||, [§]: 777,1/2 < hundreds of GeV, whereas m < (2 — 3) TeV. 

In the present paper we will vary the parameters of the minimal SUGRA scheme in wide 
ranges in order to carefully analyze the impact of the low reheating-temperature cosmology 
on relic neutralinos. The ranges we adopt are the following: 

50 GeV < mi/2 < 3 TeV , 

m <3TeV, (17) 
-3< A <+3, 
1 < tan/? < 60. 

The sign of [i is chosen to be positive, since negative values are somehow disfavoured by the 
limits on the supersymmetric contribution to the muon anomalous magnetic moment. On 
the configurations obtained by randomly scanning the above defined parameter space, we 
apply the experimental limits quoted above on Higgs and supersymmetry searches and on 



the b — > s + 7 decay (for details, see for instance Ref . p2|. p3| ). 

We finally remark that the phenomenology of relic neutralinos, in some sectors of the 
minimal SUGRA scheme, is also quite sensitive to some Standard Model parameters, such 
as the top quark mass m t , the bottom quark mass nib and the strong coupling a s || ^3|, |24| . 
For these parameters, we use here their 95% CL ranges: mf ole = (175 ±10) GeV, m^Mz) = 
(3.02 ± 0.21) GeV and a s (M z ) = 0.118 ± 0.004. 
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IV. RESULTS AND DISCUSSION 



As already discussed in Sec.|J, the calculation of the neutralino relic abundance Q x h 2 , 
both in the standard case outlined in Sec. |U A| and in the low-reheating scenario described 
in Sec. [II B| , relies on a detailed calculation of the neutralino self-annihilation cross section 
(o'a.imv). Here we consider the universal SUGRA model outlined in the previous Section, and 
calculate (cr ann t>) following the procedure given in Ref. [19|, to which we refer for details. 

A first preliminary conclusion about neutralino dark matter and models with a low- 
reheating temperature may be drawn by calculating the neutralino annihilation cross section 
at the temperature T*. The quantity (cr ann t>)* is shown in Fig.|I| as a function of m x in terms 
of a scatter plot obtained by varying the SUGRA parameters in the ranges given by Eq. 
(|T%|). The scatter plot is compared to the values of the limiting cross section given in Eq. 
CpTT|) , calculated for different values of T RH . This figure shows that the condition of non- 
relativistic non-equilibrium, given by Eq. ([11]), is verified only for values of Tr# smaller 
than about 300 MeV. For larger values of the reheating temperature, (cr arm t>)* always lies 
above the curves of the limiting cross section. This fact implies that for Trh > 300 MeV 
the neutralino always reaches thermal equilibrium during the reheating phase, and therefore 
the peculiar behaviour Q x h 2 oc (o" ann t>) is limited only to cosmological models with very low 
reheating temperatures. This conclusion is actually true for a class of supersymetric scenarios 
which is more general than the universal SUGRA model shown in Fig.[l|. For instance, we 
have explicitly verified that the same result also applies to other SUGRA schemes where 
the universality condition of Eq. (pT5|) are relaxed at the GUT scale for the Higgs sector, 
when the supersymmetric parameters are varied in the same ranges discussed in the previous 
Section. One must notice that the effective lower limit on neutralino cross section revealed 
in Fig. |l| is a consequence of the choice of the upper ranges of Eq. (|18|) adopted for the 
dimensional supersymmetric parameters. However, we remind that the intervals given in 
Eq. fll^D are representative of a typical upper bound on the supersymmetry breaking scale, 
M SUBy < a few x 10 3 GeV. This constraint may be understood on quite general grounds, 
since it derives from naturalness arguments on the stability of the Higgs potential and the 
requirement of absence of fine-tuning in the generation of the electroweak scale through the 
mechanism of radiative symmetry breaking. Notice that, as far as m\/2 is concerned, the 
upper value adopted in Eq. (|18D is already large enough to pose questions about fine-tuning 
problems |3|, |J]. 

Let us turn now to the calculation of the neutralino relic abundance. Fig.0 shows Q x h 2 
as a function of the neutralino mass m x , calculated in universal SUGRA for tan/3 = 30 
and Aq = and with mo and mj/2 varied in the ranges of Eq. fll8[) . Starting from the 
upper-left panel, we plot the results obtained for standard cosmology, and for low reheating- 
temperature cosmologies with Trh = 30, 20, 10, 5, 1 GeV. In the case of standard cosmology, 
the neutralino relic abundance turns out to be generically quite large, in excess of the upper 
bound on the total amount of matter in the Universe given in Eq. (U). This is a typical 
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feature of universal SUGRA models for values of tan /3 < 40 [g, [§ g, §, % |, 0, [T|, |T|, [L3| . 
For this schemes, the constraint on the parameter space coming from cosmology is actually 
very strong, especially in posing stringent upper limits on the neutralino mass. Fig. |2| shows 
that an upper bound of about 200 GeV is obtained on the neutralino mass when we restrict 
Q x h 2 to be less than 0.3. 

There are, however, ways out to avoid such a bound. When coannihilation is included, 
the upper bound on m x in universal SUGRA at low/intermediate tan (3 can be extended up 
to about 500 GeV || [7|, [11], |13| , even though this possibility is restricted to a very narrow 
sector of the SUGRA parameter space where the neutralino mass is almost degenerate with 
the stau mass. 

Alternatively, when a low reheating temperature is allowed, Fig.|] shows that the upper 
limit on m x coming from cosmology is removed. By lowering the reheating temperature we 
affect mainly the relic abundance for large neutralino masses, unless we lower T R h below a 
few GeV: for T R h < 5 GeV, £l x h 2 is reduced for all the allowed mass range of the SUGRA 
model. 

The behaviour of the different panels of Fig.|2| may be easily understood by comparing the 
relevant mass scales which enter in the calculation of the relic abundance in the standard and 
low reheating-temperature models, namely m x , Tp and T R h- In Fig-H we show a generic 
example of what happens at different neutralino masses for a sufficiently low reheating 
temperature. The upper panel of Fig.|3] shows the relic abundance vs. m x , calculated for 
(oarm^) fixed at the value 10~ 10 GeV -2 . The thick horizontal (black) line refers to the 
calculation for standard cosmology, while the thick decreasing (blue) line refers to T RH = 10 
GeV. The two lines are superimposed for m x < 150 GeV. The lower panel shows the value 
of the freeze-out temperature Tp as a function of the neutralino mass. The horizontal line 
indicates the values of T R h = 10 GeV. By comparing the upper an lower panels, we see 
that as long as the freeze-out occurs at a temperature which is (much) smaller than Trh, no 
difference is present for the low reheating-temperature scenario as compared to the standard 
case. On the contrary, when T F approaches T RH , the relic abundance becomes suppressed, 
and this effect becomes more pronounced as T F grows as compared to T RH . Since T F is an 
increasing function of the neutralino mass, a low-reheating scenario affects more the large 
mass sector of the theory, and this explains the peculiar behaviour of Fig. 0. 

Going back to the discussion of the features of Fig. 0, we see that in order to have 
heavy neutralinos compatible with the cosmological upper bound in Eq. ([|), the reheating 
temperature has to be lowered below about 30 GeV. When T R h falls below 20 GeV all the 
neutralino mass range is acceptable from the point of view of cosmology, although a fraction 
of the SUGRA configurations entail relic neutralinos which are a subdominant dark matter 
component, since Q x h 2 falls below the lower limit of Eq. (g). When T RH < 0.6 GeV all the 
SUGRA configurations lead to subdominant relic neutralinos. For smaller values of T R h the 
neutralino relic abundance is even more suppressed and becomes negligible for T R h in the 
MeV range. 
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We comment at this point that for larger values of tan /3 (tan j3 > 40) in SUGRA models an 
upper limit on the neutralino mass is not present even in standard cosmology. The occurrence 
of an acceptable relic abundance is generically confined to corridors in the parameter space 
where either the relic abundance is suppressed by coannihilation or the neutralino mass lies 
close to the pole of the annihilation cross section mediated by the pseudoscalar higgs field 
A. In a low reheating-temperature cosmology, these features are relaxed and conclusions 
similar to those discussed above in connection with Fig.||] are present. In particular, a scale 
of about 1 GeV as a lower limit on T R h in order to have dominant relic neutralinos is present 
for all the values of tan f3 given in Eq. (|18|) . 

New features in the parameter space are also present when the universality is not ex- 
tended to Higgs masses H, U, H, [TTH . Effects of non-universality modify the low-energy 



sector of the theory, through the RGE evolution and the conditions of radiative electroweak 
symmetry breaking: this alters the neutralino couplings and the mass spectrum of sparticles 
and induces variations of the neutralino annihilation cross section. The consequence is a 
change and an extension of the regions of parameter space which are compatible with a relic 
neutralino, also in standard cosmology. A low reheating temperature has again the effect 
of enlarging these cosmo logically relevant sectors of the parameter space. A lower limit of 
about 1 GeV on Trh in order to have dominant relic neutralinos is again recovered also in 
non-universal SUGRA models. 

Fig. H shows the same information contained in Fig. |2|, expressed in the plane mo-^i/2- 
The upper-left panel refers to standard cosmology, the other panels refer to low reheating- 
temperature cosmologies with T RH = 30, 20, 10, 5, 1 GeV. The value of tan/3 is fixed at 30 
and A = 0. The light-shaded (yellow) regions are excluded domains: the ones on the left 
side of the m -mi/ 2 plane are excluded by the experimental bound discussed in the previous 
Section or by the non-occurrence of radiative electroweak symmetry breaking, while the ones 
on the lower part of the plane do not correspond to viable LSP neutralino models. The dark- 
shaded (blue) areas correspond to the domains where the neutralino relic abundance falls 
inside the cosmologically relevant range of Eq. (fj). The hatched (red) regions correspond 
to neutralinos with a subdominant relic abundance, i.e. to Q x h 2 < 0.05. 

The first panel shows that for standard cosmology the region of the m -mi/2 parameter 
space which is allowed by cosmology is quite restricted, as discussed above. The upper 
limit on the acceptable values of the neutralino relic abundance poses severe bounds on 
both itlq and m^, of the order of a few hundreds of GeV. When coannihilation is included 
f|. |7|, [ll], [12], pL3| , the allowed regions is extended in a thin band close to the lower excluded 
area. However, as is evident from the other panels, when the reheating temperature is 
lowered the consequent suppression on Q x h 2 weakens considerably the constraints on m 
and mi/2- In particular, the regions of the SUGRA parameter space which are compatible 
with the assumption of a dominant neutralino dark matter change significantly, depending 
on the actual value of T RH . For values of T RH below about 20 GeV a large fraction of 
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the parameter space is allowed by cosmology and for Trh < 1 GeV all the supersymmetric 
parameter space is compatible with the cosmological abundance of relic neutralinos, even 
though for most of the values of the parameters the relic abundance falls below the interval 
ofEq. (|). 

As we have discussed before, the extension of the cosmologically allowed regions in the 
plane m -mi/ 2 is enlarged, in standard cosmology, when tan (3 is larger than about 40, since 
in this case coannihilation or annihilation through the A-pole is effective in reducing the 
values of Q x h 2 . In the latter case, an almost diagonal allowed band opens up || ^4|. Larger 
regions are also allowed in non-universal SUGRA models. In both low reheating 

temperature has again the effect of widely enlarging the cosmologically relevant domains 
in the plane mQ-mi/2. We obtain, also for the SUGRA models with tan/? > 40 and for 
the non-universal SUGRA schemes, that in all the supersymmetric parameter space the 
neutralino relic abundance is compatible with the range of Eq. (0) when Trh < 1 GeV. 

In the above discussion, we noticed that, for any given neutralino mass, low enough values 
of T RH would imply a suppression of fi x h 2 too strong to be compatible with the hypothesis 
of dominant neutralino dark matter. We can therefore use this argument to obtain a lower 
limit on Trh under the assumption that the neutralino represents the dominant component 
of dark matter in the Universe. This lower bound on Trh is shown as a solid line in Fig|| 
as a function of m x . To derive this limit we have varied all the SUGRA parameters in the 
intervals of Eq. ([TJJ). The result plotted in Fig. |5] shows that if we require the neutralino to 
be the dominant dark matter component, the reheating temperature in the early Universe 
cannot be lower than a value ranging from 0.6 GeV up to about 20 GeV, depending on the 
value of the neutralino mass. This bound on T RH at the GeV scale is quite interesting, since 
it is sizably stronger than the constraint given by nucleosynthesis. We have verified that the 
same result remains valid also for non-universal SUGRA models. Therefore a lower limit of 
about 1 GeV on Trh, in order to explain the dark matter content of the Universe in terms 
of relic neutralinos, is a specific feature of supergravity models. 



V. CONCLUSIONS 

In standard cosmology it is usually assumed that the temperature Trh of the Universe 
at the beginning of the radiation-dominated era is much higher than the supersymmetry 
breaking scale. Moreover, neutralinos decouple from the thermal bath after the reheating 
phase, which followed the end of inflation, has terminated. Under these assumptions, the 
allowed parameter space of SUGRA models turns out to be severely constrained by the 
requirement that the neutralino relic density does not exceed the maximal value of the 
matter density of the Universe deduced from observations. 

In this paper we have performed a quantitative study of the neutralino relic abundance 



in cosmological scenarios with a low reheating temperature [ TA \ . This is a viable possibility 
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since the only robust lower bound on the reheating temperature Trh can be set at the 
MeV scale, in order not to spoil nucleosynthesis predictions. The suppression on Q x h 2 is 
originated by the fact that neutralinos decouple from the thermal bath before the end of the 
reheating phase. In this case the neutralino number density is diluted by entropy production 
and by a higher expansion rate than in the radiation-dominated era. 

For values of Trh < 30 GeV the domains of the SUGRA parameter space which are 
compatible with dominant relic neutralinos are largely enhanced with respect to the standard 
cosmological case. These domains depend on Trh and we have shown their evolution as a 
function of the reheating temperature. For Trh < 1 GeV all the SUGRA parameter space 
becomes compatible with the bounds on the dark matter relic abundance (even though the 
neutralino relic abundance for these low values of Trh is strongly suppressed) . 

Since lower Trh imply smaller relic densities, the assumption that neutralinos provide a 
major contribution to the dark matter of the Universe implies a lower limit on Trh- This 
constraint ranges from 0.6 GeV for neutralino masses of the order of few tens of GeV, up to 
20 GeV for neutralino masses in the TeV range. This bound on T RH , subject to the request 
of explaining the dark matter content of the Universe only in terms of relic neutralinos in 
SUGRA schemes, is much stronger than the limit on Trh coming from nucleosynthesis. 
Similar conclusions occur also for non-universal SUGRA models. 
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FIG. 1: Thermal average, at the temperature ~ m-y/A, of the neutralino self-annihilation cross 
section times the relative velocity (<7 ann v)* as a function of the neutralino mass m x . The points 
denote the values of (<T ann v)* calculated in universal SUGRA with the parameters varied as in 
Eq. (|l^). The values of mt° le , rrib and a s are varied inside their 2a allowed intervals. The solid 
lines denote, for different values of the reheating temperature Tr#, the values of the limiting 
cross section of Eq. ( |Tl| ) which determines, for a cosmological model where neutralinos are always 
non-relativistic, whether neutralinos can reach thermal equilibrium during the reheating phase. 
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FIG. 2: Neutralino relic density 0, x h 2 as a function of the neutralino mass m x in universal SUGRA, 
for tan/3 = 30 and Aq = 0. The parameters rriQ and m\ii are varied according to the intervals 
of Eq. dl§|). The values of mf ole , mb and a s are fixed at their central values: mf ole = 175 GeV, 
rrib(Mz) = 3.02 GeV and a s (Mz) = 0.118. The upper-left panel shows results for standard 
cosmology. The other panels refer to different values of the reheating temperature Trh- The 
horizontal solid and dashed lines delimit the interval for the amount of non-baryonic dark matter 
in the Universe, given in Eq. (Q). 
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FIG. 3: Effect of a low reheating-temperature on the WIMP relic abundance. The upper panel 
shows ^ x /i 2 as a function of m x for standard cosmology (upper thick horizontal line) and for 
Trh = 10 GeV (thick decreasing line). The two lines are superimposed for m x < 150 GeV. The 
thin horizontal solid and dashed lines delimit the interval for the amount of non-baryonic dark 
matter in the Universe, given in Eq. (Q). The lower panel shows the freeze-out temperature Tp as 
a function of m x . The horizontal line denotes the values of the reheating temperature: Trh = 10 
GeV. 



19 




T.„ = 30 GeV 



500 1000 1500 2000 2500 3000 

m,„ (CeV) 




2500 3000 




500 1000 1500 2000 2500 3000 



500 1000 1500 2000 2500 3000 



(GeV) 



(GeV) 





FIG. 4: Cosmologically favoured regions in the plane m±/2~ m o i n universal SUGRA for tan/3 = 30 
and Aq = 0. The values of m^ ole , rrif, and a s are fixed at their central values: m^ ole 



175 

GeV, rrib(Mz) = 3.02 GeV and a s (Mz) = 0.118. The upper-left panel shows results for standard 
cosmology. The other panels refer to different values of the reheating temperature Tr#. The 
dark-shaded (blue) areas correspond to configurations where 0.05 < ^x^ 2 — 0-3. Hatched (red) 
areas denote configurations where ^ x /i 2 < 0.05. Light-shaded (yellow) regions are excluded either 
by experimental constraints or by theoretical arguments. 
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FIG. 5: Lower limit on the reheating temperature Trh as a function of the neutralino mass m x , 
obtained by requiring that the neutralino is the dominant component of dark matter in the Universe 
(i.e.: 0.05 < f^ x /i 2 < 0.3). The result refers to universal SUGRA with the parameters varied as in 
Eq. (|l~8|). The values of mf ole , and a s are varied inside their 2a allowed intervals. 
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